Isotropic diffusion of the small ribosomal subunit in Escherichia coli 
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The ribosome is one of the most important macromolecular complexes in a living cell. We 
have tracked individual 308 ribosomal subunits in exponentially growing E. coli cells using three- 
dimensional single-particle tracking, where the z-position is estimated using astigmatism. The 30S 
subunits are stochiometrically labeled with S2-mEos2 by replacing the rpsB gene in the E. coli 
chromosome by rpsB-mEos2. The spatial precision in tracking is 20 nm in xy and 70 nm in z. 
The average trajectory consists of 4 steps corresponding to 80 ms. The trajectories are excluded 
from parts of the cell, consistent with nucleoid exclusion, and display isotropic diffusion with nearly 
identical apparent diffusion coefficients (0.05 =b 0.01 ^im^ s~^) in y and z. The tracking data fits 
well to a two-state diffusion model where 46% of the molecules are diffusing at 0.02 \im^ s~^ and 
54% are diffusing at 0.14 [irn^ s~^. These states could correspond to translating 70S ribosomes and 
free 30S subunits. 



Fluorescence microscopy has revolutionized cell biol- 
ogy during the last decades. However, a fundamental 
limitation of using optical microscopy in subcellular stud- 
ies in bacteria has been the inherent resolution limita- 
tion of approximately 200 nm, imposed by diffraction. 
This makes it impossible to distinguish two objects that 
are closer than 200 nm. For fluorescence microscopy, 
diffraction-unlimited resolution is achievable by inhibit- 
ing the simultaneous fluorescence of fluorophores within a 
diffraction-limited area, i.e. by spacing the observations 
in time instead of space. One group of techniques are 
capable of achieving this by stochastically allowing fluo- 
rescence from a small subset of fluorophores and are com- 
monly referred to as (F)PALM or STORM [1^^. These 
techniques rely on the minute probability of activating 
multiple fluorophores within a diffraction-limited area. 
To find the spatial position centroid calculation has to be 
performed and this is most commonly achieved by fitting 
a two-dimensional Gaussian function to the point spread 
function (PSF). Since these stochastic methods are inher- 
ently easy to integrate in a normal wide-field fluorescence 
microscope with basically only a strict requirement of flu- 
orophores that can be converted or switched on and off, 
they are widely used in biology. Since the development 
of photoconvertible fluorescent proteins, biological appli- 
cations for super-resolution imaging have started to ap- 
pear and are becoming increasingly popular due to their 
non-invasive nature [3 . Lately these techniques have ex- 
panded from being inherently two-dimensional to yield- 
ing position data in all three dimensions. This is most 
commonly achieved by either introducing an optical 'de- 
fect', astigmatism, in form of a cylindrical lens in front 
of the camera (astigmatism approach) [4 , or by splitting 
the fluorescence and detecting it on two different parts of 
the camera with one having a longer light path, leading 
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to images of different focal planes (biplane approach) [5] . 
More exotic versions based on for example interferome- 
try [6 , 4Pi setups [7 and complex PSF shaping setups 
relying on e.g. adaptive optics [8 are also available. 

Single-particle tracking (SPT) with nanometer preci- 
sion (~10 nm) for single molecules was initially demon- 
strated by tracking individual lipids in a lipid membrane 
[9 . The particle tracking was extended to the third di- 
mension by Kao and Verkman who introduced three- 
dimensional (3D) tracking using astigmatism An- 
other step forward was taken by Manley et al. and Niu et 
al. who combined SPT with stochastic super-resolution 
techniques by using photoconvertible proteins [TTJ [T2] . 
This made it possible to work with abundant molecules 
by converting a few molecules at a time. When the 
molecules bleach more molecules can be converted and 
tracked in the same cell until a sufficient number of tra- 
jectories have been acquired. Here we have combined 
tracking of photoconvertible proteins with 3D localiza- 
tion to study the diffusion of the small ribosomal subunit 
labeled with the photoconvertible protein mEos2 fused to 
the ribosomal protein S2. 



RESULTS AND DISCUSSION 

We obtained 1222 S2-mEos2 trajectories from one Es- 
cherichia coli {E. coli) cell (see Figure [T]). The trajecto- 
ries form a cylinder with hemispherical end caps, which 
is the typical shape of an E. coli cell. Figure |2] shows 
two individual S2-mEos2 trajectories which are confined 
to regions with radii of approximately 100 and 150 nm. 
This confinement is most likely due to the tethering of 
the 30S subunit to mRNAs [13]. 

To investigate whether the diffusion of the small sub- 
unit is isotropic or not, we calculated the y and z mean 
square displacements (MSDs) of the trajectories (see Fig- 
ure H). The MSD curves plateau after 200 ms and the 
X and y MSDs are nearly identical, indicating that cellu- 
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FIG. 1: Overlay of 1222 S2-mEos2 trajectories obtained from 
one E. coll cell. The trajectories form a cylinder with hemi- 
spherical end caps, which is the typical shape of an E. coli 
cell. 
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FIG. 2: Two S2-mEos2 trajectories obtained from one E. coli 
cell. The trajectories are confined to regions with radii of 
approximately 100 and 150 nm. 



lar confinement has a very small effect on S2-mEos2 dif- 
fusion at this timescale. The z MSD curve has a higher 
offset than the x and y MSD curves but is otherwise 
very similar to these curves. The higher offset can be 
explained by the higher uncertainties of the 2;-positions. 
The X, y and z apparent diffusion coefficients obtained 
by fitting lines to MSDs up to 100 ms are 0.04, 0.04 and 
0.06 |LLm^ s~^, respectively. Unlike the MSDs themselves, 
these coefficients are very similar since they are largely 
unaffected by the localization errors which mainly affect 
the offsets. 

To further analyze the trajectories, we fitted cumula- 



FIG. 3: Mean square displacements (MSDs) of the S2-mEos2 
trajectories shown in Figure [l] Error bars represent standard 
error of the means. All MSD curves plateau after 200 ms. The 
X and y MSDs are nearly identical while the z MSD curve has 
a higher offset than the x and y MSD curves but is otherwise 
very similar to these curves. The x, y and z apparent diffusion 
coefficients obtained by fitting lines to the first 5 MSDs are 
0.04, 0.04 and 0.06 [im^ s~\ respectively. 



tive distribution functions (CDFs) corresponding to one-, 
two-, and three-state diffusion models to the experimen- 
tal CDF of 20- ms displacements (see Figure [4|. It is clear 
that a one-state model is not sufficient in describing the 
data. The sum of the squared errors decrease two orders 
of magnitude when the number of states is increased from 
one to two but only a factor of two when the number of 
states is increased from two to three. This indicates that 
a model with two diffusion states is sufficient to explain 
the data. The apparent diffusion coefficient obtained 
from fitting the experimental CDF with a CDF corre- 
sponding to a one-state diffusion model is 0.07 iitm^ s~^. 
This apparent diffusion coefficient is higher than the 
apparent diffusion coefficients obtained from the initial 
slopes of the MSD curves, 0.05 ± 0.01 |Ltm^ s~^, since it 
also includes the localization errors. 

The apparent diffusion coefficients obtained from fit- 
ting the experimental CDF with a CDF correspond- 
ing to two diffusion states are 0.14 jitm^ s~^ for the 
fast state and 0.03 iitm^ s~^ for the slow state. When 
we take the localization errors into account we obtain 
Di = 0.14 M^m^ s-^ and D2 = 0.02 ^tm^ s-^ The occu- 
pancies of the two states are 54% for the fast state and 
46% for the slow state. The diffusion states probably cor- 
respond to free 30S ribosomal subunits and translating 
70S ribosomes. If this is the case, it implies that 54% of 
the small subunits are involved in translation under our 
experimental conditions. It should be noted that free S2- 
mEos2 proteins in the cell diffuse too fast to be localized 
with 20-ms exposures. 

The 3D tracking reveals that the ribosomes move 
isotropically in y and z. This is advantageous since 
it implies that most information about ribosome move- 
ment can be extracted from diffusion along the bacterial 
long axis (x) alone. The advantage of studying diffusion 



3 




50 100 150 200 



xy displacement (nm) 

FIG. 4: Cumulative distribution function (CDF) of the xy 
displacements over 20 ms, calculated from the S2-mEos2 tra- 
jectories shown in Figure [l] The CDF is fitted with two CDF 
models of Brownian motion corresponding to one and two 
diffusion states. The sums of the squared errors are 0.4 and 
0.004, respectively. 

along the long axis is that the diffusion in y and z will be 
spatially restricted by the membrane and possibly also 
the nucleoid. Furthermore, the confinement in y will de- 
pend on the 2;-position, which will always be the most 
inaccurate in a 3D imaging system based on astigmatism. 

To check if the small subunit is excluded from parts 
of the cell, we constructed a super-resolution image from 
all independent S2-mEos2 localizations in a 200-nm thick 
2;-section of the cell (see Figure [5|. S2-mEos2 is excluded 
from certain parts of the cell, consistent with nucleoid ex- 
clusion. This effect has previously been observed by elec- 
tron microscopy of fixed cells [14] and super-resolution 
imaging of Sra-mEos2 [15]. 

An ideal use of single-molecule tracking data is the 
determination of rates by which individual molecules go 
between different complexes with different diffusion prop- 
erties. In our case, such states could for example corre- 
spond to free small subunits and translating 70S ribo- 
somes. However, in our case, the average trajectory con- 
sists of 4 steps, corresponding to 80 ms, which implies 
that we would expect that less than 1% of the tracked ri- 
bosomes complete a round of translation within a trajec- 
tory, assuming that it would take at least 10 s to trans- 
late a typical open reading frame [16 . The number of 
translation initiation events will be just as unlikely in a 
steady-state situation. Therefore, in order to study the 
dynamics of translation initiation, elongation, termina- 
tion and ribosome recycling using single-molecule track- 
ing, the trajectories would have to span longer times. 
One way of achieving this is by spacing out the imag- 
ing events such that the molecules can be monitored for 
longer times at the same total exposures, since photo- 
bleaching is often the limiting factor. The disadvantage 
of spacing out the observation time points is that it is 
less certain that the same molecule is being tracked, i.e. 
the molecule that is being imaged can photobleach in 
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FIG. 5: Bright-field image of an E. coli cell expressing S2- 
mEos2, a super-resolution image constructed from all in- 
dependent S2-mEos2 localizations in the cell and a super- 
resolution image constructed from all independent S2-mEos2 
localizations in a 200-nm z section of the cell. Note that 
S2-mEos2 is excluded from certain parts of the cell, which is 
consistent with nucleoid exclusion. 



one frame while another molecule spontaneously pho- 
toconverts before the next imaging event. The rate of 
spontaneous photoconversion of an mEos2 molecule is 
not known but can be considered to be very low. How- 
ever, since there are so many ribosomes per cell (^25000) 
[16 , the overall risk of spontaneous conversion under a 
time period corresponding to a translation cycle might 
be significant. 

To conclude, the 30S subunits are excluded from parts 
of the cell, consistent with nucleoid exclusion, the aver- 
age apparent diffusion coefficient is 0.05 iitm^ s~^, and 
it appears that S2-mEos2 diffuses in at least two differ- 
ent complexes that are likely to be free 30S subunits and 
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translating ribosomes. 

METHODS 
Optical setup 

The optical setup included a 405-nm photoconver- 
sion laser (Radius 405-50, Coherent), a 555-nm exci- 
tation laser (GCL-150-555, CrystaLaser), a cylindrical 
lens (LJ1516RM-A, Thorlabs), an EMCCD camera (DU- 
897E-CS0-#BV, Andor Technology), an emission filter 
(HQ605/75m, Chroma Technology), an excitation fil- 
ter (Z550/20x, Chroma Technology), an image split- 
ter (Optosplit II, Cairn Research), an inverted micro- 
scope (Eclipse Ti-E, Nikon Instruments), two dichroic 
mirrors (Z405RDC and Z555RDC, Chroma Technology), 
two shutter drivers (T132 and VMM-D3, Vincent Asso- 
ciates) and two shutters (LS6T2 and LS6ZM2, Vincent 
Associates). MetaMorph 7.7.6.0 (Molecular Devices) was 
used to control the LS6T2 shutter, the camera and the 
microscope. 



Strain construction 

A strain expressing S2-mEos2 with a glycine linker 
was constructed by lambda Red recombination [18 of 
a rpsB-mEos2 DNA fragment in BW25993 containing 
the pKD46 plasmid. The construct was transferred to 
a clean BW25993 background by PI phage transduction 
and the chloramphenicol resistance cassette was removed 
by expressing FLP from pCP20. 

Sample preparation and fluorescence imaging 

The cells were grown overnight in M9 minimal medium 
supplemented with 0.4% glucose and RPMI 1640 amino 
acids (R7131, Sigma-Aldrich) at 37°C, diluted 1:1000 in 
the morning and grown in the same manner to an optical 
density at 600 nm of 0.2 and placed on a 2.5% agarose 
pad containing fresh medium. The cells were incubated 
for 30 min at room temperature (21° C) and subsequently 
imaged for 2 min at 50 Hz. The laser power density of 
the excitation laser (555-nm) was 1 kW cm~^. 



Calibration of photon numbers 

The baseline offset and the standard deviation of the 
camera noise were determined by acquiring 200 images 
with 256 X 256 pixels at 50 Hz with a closed shutter, and 
calculating and fitting a Gaussian function to the ex- 
perimental probability mass function of the intensities. 
The multiplication factor was determined by acquiring 
200 images with 256 x 256 evenly illuminated pixels at 
50 Hz with a white paper placed in front of the objec- 
tive, flattening them to correct for any unevenness in 
the illumination, and calculating and fitting a theoretical 
probability mass function of intensities to the experimen- 
tal probability mass function of the intensities [17 . The 
baseline offset, the multiplication factor and the quantum 
efficiency of the camera were used to convert the image 
intensities from counts to number of photons. 



Image analysis 

PSFs of S2-mEos2 molecules were identified with a 
wavelet segmentation algorithm [19] and fitted with 
two-dimensional Gaussian functions. For each detected 
molecule the localization uncertainty was determined by 
the fitting error while the number of photons was deter- 
mined by calculating the volume under the fitted Gaus- 
sian surface. The 2;-position was determined by minimiz- 
ing the norm given in |4]. Points were accepted if the 
fit was successful, the norm was less or equal to 5, the 
number of background photons was greater than zero, 
the number of photons within the Gaussian fit was at 
least 100, the standard uncertainty of the x-position was 
less or equal to 60 nm, the standard uncertainty of the 
^-position was less or equal to 40 nm and the 2;-position 
was between -200 and 700 nm. 



Calibration of 2;-position 

Fluorescent 40-nm beads (F8794, Molecular Probes) 
were attached to a Poly-L-lysine-treated coverslip which 
was placed on a 2.5% agarose pad (SeaPlaque GTG 
Agarose, Lonza) containing M9 minimal medium sup- 
plemented with 0.4% glucose. Individual beads were im- 
aged at 2;-positions ranging from -700 nm to 700 nm with 
50-nm increments. Point spread functions (PSFs) were 
detected with a wavelet segmentation algorithm and fit- 
ted with two-dimensional Gaussian functions. Calibra- 
tion curves were obtained by fitting cubic functions to 
the X and y widths of six beads. 



Trajectory construction 

The trajectories were constructed such that they only 
included points that could not be part of multiple trajec- 
tories. The maximum allowed x, y and z displacements 
were 200 nm, 200 nm and 300 nm, respectively. 

ACKNOWLEDGEMENTS 

We thank Prune Leroy for helpful discussions and com- 
ments on the manuscript. This work was supported 
by Vetenskapsradet, the Goran Gustafsson Stiftelse, the 



5 



Knut and Alice Wallenberg Foundation and the Foun- 
dation for Strategic Research. The authors declare that 
they have no competing financial interests. 

AUTHOR CONTRIBUTIONS 

J.E. designed the research. A.S. constructed the bac- 
terial strain. A.S. and F.P. built the optical setup. A.S. 
performed the experiments, analyzed the images and did 
the MSD analysis. F.P. did the CDF analysis. A.S., F.P. 
and J.E. discussed the analysis methods and the results, 
and wrote the paper. 



[1] M. J. Rust, M. Bates, and X. Zhuang. Sub-diffraction- 
limit imaging by stochastic optical reconstruction mi- 
croscopy (STORM). Nat Methods, 3(10):793-796, 2006. 

[2] S. T. Hess, T. P. Girirajan, and M. D. Mason. Ultra- 
high resolution imaging by fluorescence photoactivation 
localization microscopy. Biophys. J., 91(ll):4258-4272, 
2006. 

[3] Bo Huang, Hazen Babcock, and Xiaowei Zhuang. Break- 
ing the diffraction barrier: super-resolution imaging of 
cells. Cell, 143(7): 1047-1058, 2010. 

[4] B. Huang, W. Wang, M. Bates, and X. Zhuang. Three- 
dimensional super-resolution imaging by stochastic opti- 
cal reconstruction microscopy. Science, 319(5864) :8 10- 
813, 2008. 

[5] M. F. Juette, T. J. Gould, M. D. Lessard, M. J. 
Mlodzianoski, B. S. Nagpure, B. T. Bennett, S. T. Hess, 
and J. Bewersdorf. Three-dimensional sub-100 nm res- 
olution fluorescence microscopy of thick samples. Nat 
Methods, 5(6):527-529, 2008. 

[6] G. Shtengel, J. A. Galbraith, C. G. Galbraith, 
J. Lippincott- Schwartz, J. M. Gillette, S. Manley, 
R. Sougrat, C. M. Waterman, P. Kanchanawong, M. W. 
Davidson, R. D. Fetter, and H. F. Hess. Interferomet- 
ric fluorescent super-resolution microscopy resolves 3D 
cellular ultrastructure. Proc. Natl. Acad. Sci. U. S. A., 
106(9):3125-3130, 2009. 

[7] M. Nagorni and S. W. Hell. 4Pi-confocal microscopy pro- 
vides three-dimensional images of the microtubule net- 
work with 100- to 150-nm resolution. J. struct. Biol., 
123(3):236-247, 1998. 



[8] S. R. Pavani, M. A. Thompson, J. S. Biteen, S. J. 
Lord, N. Liu, R. J. Twieg, R. Piestun, and W. E. Mo- 
erner. Three-dimensional, single-molecule fluorescence 
imaging beyond the diffraction limit by using a double- 
helix point spread function. Proc. Natl. Acad. Sci. U. S. 
A., 106(9) :2995-2999, 2009. 

[9] T. Schmidt, G. J. Schtz, W. Baumgartner, H. J. Gruber, 
and H. Schindler. Imaging of single molecule diffusion. 
Proc. Natl. Acad. Sci. U. S. A., 93(7):2926-2929, 1996. 
[10] H. P. Kao and A. S. Verkman. Tracking of single fluores- 
cent particles in three dimensions: use of cylindrical op- 
tics to encode particle position. Biophys. J., 67(3): 1291- 
1300, 1994. 

[11] S. Manley J. M. Gillette, G. H. Patterson, H. Shroff, 
H. F. Hess, E. Betzig, and J. Lippincott- Schwartz. High- 
density mapping of single- molecule trajectories with 
photoactivated localization microscopy Nat. Methods, 
5(2):155, 2008. 

[12] L. Niu and J. Yu. Investigating intracellular dynamics of 
FtsZ cytoskeleton with photoactivation single-molecule 
tracking. Biophys. J., 95(4):2009-2016, 2008. 

[13] B. P. EngUsh, V. HauryUuk, A. Sanamrad, S. Tankov, 
N. H. Dekker, and J. Elf. Single-molecule investigations 
of the stringent response machinery in living bacterial 
cells. Proc. Natl. Acad. Sci. U. S. A., 108(31):E365, 2011. 

[14] C Robinow and E Kellenberger. The bacterial nucleoid 
revisited. Microbiol. Rev., 58(2) :21 1-232, 1994. 

[15] W. Wang, G-W. W. Li, C. Chen, X. S. Xie, 
and X. Zhuang. Chromosome organization by a 
nucleoid- associated protein in live bacteria. Science, 
333(6048): 1445-1449, 2011. 

[16] H. Bremer and P. Dennis. Modulation of chemi- 
cal composition and other parameters by growth rate. 
In Fredrick C. Neidhardt, editor, Echerichia coli and 
Salmonella: Cellular and Molecular Biology. ASM press, 
Washington, 1996. 

[17] Maximilian H. Ulbrich and Ehud Y. Isacoff. Subunit 
counting in membrane-bound proteins. Nat. Methods, 
4:319-321, 2007. 

[18] K. A. Datsenko and B. L. Wanner. One-step inacti- 
vation of chromosomal genes in Escherichia coli K-12 
using PGR products. Proc. Natl. Acad. Sci. U. S. A., 
97(12) :6640-6645, 2000. 

[19] I. Izeddin, J. Boulanger, V. Racine, C.G. Specht, 
A. Kechkar, D. Nair, A. Triller, D. Choquet, M. Dahan, 
and J.B. Sibarita. Wavelet analysis for single molecule 
localization microscopy Opt. Express, 20(3) :208 1-2095, 
2012. 



